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SUMMAEY 


A one-dimensional analysis of ram- jet buzz is presented. The par- 
ticular mechanism leading to buzz, which is treated herein, is the am- 
plification of acoustic waves in the ram-jet combustion chamber due to 
successive reflections from the inlet and exit sections. 

It is assimed that the combustion chamber is of constant area and 
is long con^jared i/ith the length of the inlet. The inlet of such a con- 
figuration operates quasi steadily, or nearly so, during unsteady opera- 
tion. It is assiimed that the exit operates quasi steadily at all times. 
The configuration is shown to be unstable -vdien the real part of the 
acoustic ia^jedance of the inlet is greater than a term of the order of 
the combustion-chamber Mach number. For quasi-steady operation, the im- 
pedance of an inlet is proportional to the slope of its characteristic 
c\irve. Increasing the combustion-chamber Mach number or decreasing the 
slope of the inlet characteristic curve during subcrltlcal operation 
will tend to increase the range of stable operation. 

Cortputations indicate that burning increases the stable operating 
range of a given configuration. The computations assume a fixed planar 
flame front and a cons taint heat release per unit mass. 

The requirement of a relatively long combustion chamber is relaxed 
in the appendix that treats the case of isentropic flow in the combus- 
tion chamber. 

The applicability of the results for configurations ^herein the 
inlet is replaced by a compressor is mentioned briefly in CONCLUDING 
EEMAEKS. 


INTRODUCTION 

Self- sustained oscillations frequently occur in systems throiagh 
v/hich a fliild is flowing. The oscillations are Induced by that com- 
ponent (or components) of the system that tends to amplify incident 
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pressure waves. The surging of a system containing a centrifugal- or 
axial- flow compressor (e.g., refs. 1 to 4) is an exan^le of a self- 
sustained oscillation. Inlets designed for supersonic jet engines also 
have heen observed to induce oscillations (e.g., refs. 5 to 11 ) vdilch 
are usually referred to as "buzz." The origin of buzz in ram- jet en- 
gines is the subject of the present report. 

Consider a nose inlet and an exit section connected by a combustion 
chamber of constant area, as indicated schematically in figure 1. Such 
a configuration has the essential features of a ram- jet engine. If the 
exit section is wide open, the normal shock is downstream of the inlet 
cowl lip, and the engine is said to be operating supercritical] y (fig. 
1(a)). Decreasing the exit area tends to move the shock forward, in- 
creasing the pressure in the combtistion chamber but leaving the mass 
flow unaffected. With further decreases in exit area, the shock even- 
tually moves forward of the inlet cowl lip and the inlet is said to be 
operating subcriticelly (fig. l(b)). During subcrltical operation, the 
mass flow decreases with decreases in exit area. Eventually a mass flow 
is reached below which operation is unstable, this resulting in buzz. A 
typical plot of combustion-chamber pressure against mass flow is shown 
in figure l(c). The slope of this curve is generally positive at the 
mass flow corresponding to the transition between stable and unstable 
operation. 

The origin of b\xzz may be attributed to either of two types of phe- 
nomenon. These are distinguished, herein, as follows: 

(1) Linear instability origin: Consider a weak pressure wave to be 

propagating upstream in the ccnibustion chamber. The wave will be re- 
flected by the inlet section and will travel do-vnistream Tdiere it will be 
reflected by the exit section. If the successive reflections tend to 
increase the an^lltude of the wave, the system is unstable and oscilla- 
tions will result. If the magnitu^s of the corresponding incident and 
reflected waves at the inlet are linearly related, then the oscillation 
is termed herein as having a linear instability origin. The essential 
feature of this type of instability origin is that the ratio of the mag- 
nitudes of the incident and reflected waves at the inlet be constant, 
the constant depending on the nominal stea(3y flow through the ram jet. 
This Instability can be analyzed by the method of small pearburbations . 
The entire ram- jet configuration must be considered. 

( 2 ) Nonlinear instability origin: In seme cases, a weak incident 

wave may create large-scale aerodynamic changes within the inlet (e.g., 
large-scale separation, shock motion, etc.) such that the magnitude of 
the reflected wave is not directly related to the magnitude of the in- 
cident wave. Oscillations induced by such inlet behavior are termed 
herein as having a nonlinear instability origin. The weak incident 
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wave may be thought of as "triggering" the buzz cycle. The occurrence 
of this type of Instability is expected to depend primarily on the aer- 
odynamic behavior of the inlet without particular regard for the remain- 
der of the configuration. 

It is desirable to determine criteria defining the stable operating 
range of ram- jet configurations. In reference 7 it was observed that a 
series of nose inlets buzzed at that mass flow where a vortex sheet 
(originating at the upstream shock configuration) just enters the cowl 
lip. A qualitative analysis of the flow for this condition suggested 
that the vortex sheet causes separation on the inner surface of the 
cowling. This separation apparently triggered the btizz cycle. However, 
reference 11 presents data which indicate that the entrance of the vor- 
tex sheet into an inlet cowling need not necessarily result in buzz. In 
reference 11, it is suggested that the occurrence of buzz in the experi- 
mental ram-jet configuration test reported therein mi^t be associated 
with separation on the inlet centerbody. The role of separation and 
boundary-layer build-tp as mechanisms leading to buzz is also discussed 
in reference 8. Both references 8 and 11 present extensive experimental 
data and discussions of the nature of the buzz cycle. 

Several papers, such as references 6 and 9, have Investigated ana- 
lytically the stability of ram- jet configurations having a linear insta- 
bility origin. Eeference 6 likens the ram- jet configuration to a 
Helmholtz resonator having a mean through flow. The system beccmes un- 
stable when the slope of the characteristic curve (pressure against mass 
flow) of the inlet becomes sufficiently positive. Both burning and non- 
burning cases -vrere considered. For the nonbuming case, conditions in 
the combustion chamber eire assumed uniform at any instant, the wave 
structure in the combustion chamber thus being neglected. This assump- 
tion is valid provided that the combustion chamber is very short com- 
pared with the wavelength of the oscillations. A similar restriction 
applies to the analysis of the burning case. When the cambiistion cham- 
ber is not short compared ^/ith the wavelength of the oscillations, it is 
necessary to take into accoimt the existence of waves in the ccmbxistion 
chamber. The stability of a relatively simple configuration of this 
type is analyzed in reference 9, Tdiich considers a configuration con- 
sisting of an open-nose inlet, a relatively long constant-area combus- 
tion chamber, and a choked exit. The stability of the system is ana- 
lyzed using transfer functions (ref. 12). (ihe transfer-function con- 
cept provides an extension of classical acoustic-impedance techniques 
to systems having a mean through flow.) Eeference 9 reports that the 
configuration became unstable at flight Mach numbers above a certain 
critical value, the critical value depending on the Mach number of the 
mean flow through the combustion chamber. This theoretical trend \m.s 
confiimed by experiment. 
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A stability analysis, similar to that of reference 9 is presented 
herein. It is assumed that the ram jet has a constant-area combustion 
c hamb er and a linear instability origin. However, the inlet is not 
specialized to be an open-nose inlet. The flow is considered to be 
one-dimensional. (Some justification for this assun^Jtion can be de- 
duced from ref. 10, which correlates the wave structure in an experi- 
mental ram- jet combustion chamber with calculations based on one- 
dimensional flow theory. ) Both burning and nonburning cases are 
treated. 


ANALYSIS 

Flow through a tube of constant area is first considered. The 
fluid is assumed to be sli^tly perturbed from the condition of uniform 
through flow, and the corresponding equations of motion are noted. Pos- 
sible standing waves are indicated and the stability of flow throu^ a 
simple ram-jet configuration is then investigated. 


Basic Equations 

Consider a tube of constant area with a through flow of velocity 
u, pressure p, etc. Assume the flow to be slightly perturbed and 
represent the perturbations by dp. An, etc. The equations of motion, 
neglecting viscosity, heat conduction, and second-order terms, are 

Momentum; 


Continuity: 


Energy; 


State: 




dp Ap AT 
p “ p " T “ 


(la) 


(lb) 


(ic) 


(Id) 
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(Symbols are defined in appendix A.) Hie entropy perturbation is re- 
lated to the pressure, ten 5 )erature , and density perturbations by the 
relations 


^ = f -rf =-(r-i)f ^rf (2) 

The solution of equations (l) can be expressed in the form 

f = f[x - (a + u)t] + g[x + (a - u)t] (3a) 

r ^ = f [x - (a + u)t] - g[x + (a - u)t] (3b) 

r — = f[x - (a + u)t] + g[x + (a - u)tj - h(x - ut) (3c) 

^ = h(x - ut) . (3d) 

The functions f and g represent waves moving iq)stream and downstream, 
respectively, irf.th the speed of sound relative to the steady velocity u. 
The function h represents an entropy i^ave that is convected by the 
steady velocity u. 

It is now assumed that the perturbations are sinusoidal. (This as- 
sumption is not unduly restrictive since a nonsinusoidal oscillation may 
be considered as a superposition of sinusoidal waves . ) With the use of 
conplex notation, equations (3) may be expressed as 
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(the real part of the ri^t-hand side being imderstood) -where I', G, 
and k = cd + iA are coniplex constants- Equations (4) represent si- 
nusoidal -wa-ves with frequency to and angolification A. The mmerical 
•values of cd and A depend on the boundary conditions at the inlet 
Xj^ and the exit Xg of the tube. If A < 0, -fche -waves are damped. 

If A > Oj tne -waves arplify and the system is unstable. If A = 0, 
the -wa-ves are neither amplified nor damped and the system is in a 
state of neutral equilibrium. The case A = 0 is of particular in- 
terest since it generally defines a limit of stable operation. 


Examination of equations (4) indicates that the dimensionless com- 
plex ratios etc . , are Independent of time . This is a 


consequence of having chosen k 
in equations (4) . The ratios ^ = 0 - i-^ = 


to be the same for -the various waves 
» ^ and C 


rp Au 


Ct> Au 


will 


be used in later developments. These ratios may be termed acoustic im- 
pedance and entropy impedance, respectively, by analogy with classical 
acoustic theory (ref. 13) -1 The quantl-ty 9 is termed the acoustic 
resis-tance while ijr is -termed the acoustic reactance. Similar complex 
ratios -were used in reference 9- 


Erom equations ( 4 ) and the definitions of ^ and g. 



If 5 and 5 are known at some point Xj, then the coniplex cons-tants 
G/P and H/F can be e^qiressed in terms of the kno-wn quantities 


^In classical acoustic theory there is no mean -fchrou^ flow so that 
5 does not appear as a parameter. The parameter ^ is referred to as 
the specific acoustic impedance in ref . 13 . 
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and Equations (s) can then he written in the forms 



Thus, if ^ and C are known at one point in a tuhe, the values of ^ 
and ? everywhere in the tuhe are defined hy equations (6). If the 
boundary conations at the tube inlet and exit are specified, then k 
is so chosen that equations (s) satisfy the given boundary conditions. 

Equations (4) to (6) can also be applied to a configuration wherein 
the undisturbed flow consists of a series of uniform flows separated by 
interfaces ( such as planar flame fronts or shock waves ) . For each type 
of Interface, 'there is a corresponding set of bouadary conditions that 
must be satisfied. The boundary conditions that must be satisfied 
across a flame front are derived in appendix B. Note that each section 
of uniform flow will, in general, have different values for the con- 
stants G/F, H/F, M, and a. The, constant k, however, has the same 
value throu^out the system. 
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Ram- Jet Configuration 

The configuration investigated is of the form indicated in sketch 
(a). A combustion chamber of length L = is assumed to be of 



(a) Configuration investigated. 


constant area and to have a planar flame front at x^. At its upstream 

end, the combustion chamber is connected to an arbitrary siQjersonic in- 
let. At its do'imstream end, the combustion chamber is connected to a 
short exit section 'which is choked. The length 2e (distance from Xg 
to the throat of the exit section) is assumed very small conpaired iriLth 
L so that ieA «<!• 

The characteristic length of the inlet 2j_ -will later be assumed 
smal 1 compared 'with L (i.e., Z^/l « l) so -that the inlet may be con- 
sidered as operating quasi steadily, or nearly so, during unsteady op- 
eration. 2 The restriction Zi/L « 1 may not be applicable for prac- 
tical ram- jet configurations. This restriction is made beca\zse it sim- 
plifies the problem of determining the appropriate boundary conditions 
at x^ . The reason for choosing Zg/L f:« 1 is that the exit will at 

a1 1 times be operating quasi steadily. 


h- 

CD 

to 


^Quasi-steady operation Implies -that at each Instant, the flow 
throu^ a conponent is as though the operation "was steady. Nearly 
quasi-steady operation inplies a small departure from the quasi- 
steady condition. A component of an acoustical system operates quasi 
steadily if coZ/a « 1, -where Z is a characteristic length. It 
■will be shown that coL/a = 0(l) for the fundamen-tal mode of the ram- 
jet configurations considered. Therefore, Zi/L << 1 inplies that 
the inlet will operate quasi steadily, or nearly so, during unsteady 
operation. If Zg/L «< 1, the exit -will operate quasi steadily at 

1 times . 
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Conditions in the ’bumed-gas region {xf < x < Xg) axe indicated by 

the subscript b, ■wbile conditions in the unbumed gas region 

(xjL < X < Xf) are indicated by the subscript u. lihen the nonbuming 

case is considered, no special subscripts are used for those regions. 


Boundary Conditions 

Exit section . - The boundary condition at Xg is derived in ap- 
pendix C and may be expressed as 

^ " M(r - l)]e " ° 

Equation (?), \diich is valid provided that the exit section is operating 
quasi steadily, is assumed to apply throu^out this investigation. 

Combustion chamber . - The boundary conditions across the flame 
front are given by equations (B6a) and (B6b) of appendix B. With the 
heat release per un i t mass assumed constant (AQ =0), these equations 
become 


idiere 



(8a) 

(8b) 


The subscript f indicates that the quantities are evaluated at x^, 
the flame position. 


Inlet section . - If oiZ^/a is sufficiently small, the flow thro^Jgh 
the inlet during unsteady operation may be considered as quasi steady. 
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Then and can he expressed in terms of the slope of the charac- 

teristic curve of the inlet under steady operating conditions. If 
p = p(u) designates the combustion-chamher static pressure as a function 
of comhustion-chamher velocity during steady operation, then 



rvj± <3u 


(9a) 


?i = - (Ti - l)(^i + %) (9h) 

Equation (9h) is obtained from equation (2) hy noting that, for flows 
with constant stagnation teraperaturO, ^ — i- and 

— = , SO that - (r - 1)M. Both t. and , as 

a ^ r - 1 j^2 An T ^ ' '■i 

2 

given hy equations (9), are purely real. 

For subtly larger values of aiZ^/a (i.e., .nearly quasi-steady op- 
eration), the inlet impedance can he expressed as 


= [(^), ^ (^ 0 ) 


^diere 80^ and 8ilr^ are small and represent the departure of 0^ and 

i[r^ from the quasi-steady values of *^^du^ zero, respectively. 

A similar expression, representing a smal 1 departure from equation (9h), 
can he written for during nearly quasi-steady operation. 

A theoretical derivation for the values of and corre- 

sponding to a given inlet appears to he very difficult (even for the 
cases of quasi-steady and nearly quasi-steady flows). So far, only 
the open-nose inlet has heen treated analytically (ref. 9) . No such 
analytical development will he attempted herein. 


Nonhuming Case 

The configuration indicated in sketch (a) is assumed to he oper- 
ating without humlng. Equation (?) is assumed to apply at all times. 
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oi 

O) 

-j 

-4 


Ejqpresslng and ^ in tenna of an.d (through eqs. 

(e)) and suhstltuting Into equation. (?) yields 


-2ikL 

a(l-M2) 


1 + ^i 

1 - ^i 


M 


l+nzJ:M 


1 - 


M 


tkL 

aM(HM) 


1 - M 1 


(H) 


If and are specified, the corresponding value of k can he 

determined from equation (u) - This solution for k defines the fre- 
quency and amplification of possible standing -waves -wl-fchin the -tube. 

If A > 0, the -waves amplify and the system is Inherently uns-table. 

Another -vie-wpolnt is to -write equation (U) in -the form 

r - 1 


M + 1 cot 


= 


0 


kL 


a(l - m2). 


1 + i 


M 


2 sin 


kL 


a(l - M^) 


r 


ikL ^ 

|(r - i)cos 

kl 

K± aM(l-M2)l 

a(l - m2) 

J 


( 12 ) 


-vdiich defines the acoustic impedance at X4 in terms of k and y— , 

^1 

for a -tube of length L ha-ving a -through flo-w- of Mach number M and 
a short choked exit. If an inlet -with a known acoustic impedance is 
added at x^, then k must be determined so that the right-hand side 

of equation (l2) matches the specified Equations (u) and (l2) 

do not require l±/l‘ 1- 

It is not generally possible to obtain an explicit expression for 
k, in teams of and ^ particular solution of equations (ll) 

and (l2), applicable for a class of nonbumlng ram- jet configurations 
operating near the neutral stabill-fcy point, will now be presen-ted. (A 
solution of eqs. (ll) and (l2) for the case of Isentroplc flow, is pre- 
sented in appendix D . ) 

Consider a ram-jet configuration wi-th the properties 

1^/l « 1 


(13a) 
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A « 0 (l3b) 

M « 1 (13c) 

By virtue of equation (l3a), the inlet td.ll operate quasi steadily, or 
nearly so, duTi ng unsteady opeiration. Equation (l3b) signifies that 
the ram jet is operating near the point of inherent instahility. Under 
these conditions, it tdll he shown that 0^^ = 0 (m)j it then follot/s that 

= 0(m) and = 0(M).^ Thus, assume, first, that 

= 0(M); = 0(M) (14) 

This anticipates the result d± = 0(m). By using equations (l3) and 
(14), equations (ll) and (l2) become 


e ^ - [l + - (r - l)M + 0 (m 2)J (15a) 

M + x) + ^[j + 2N) - + 0 (m 2) (15b) 

E^ressing as ^®ij ” where 60 and 6ilrj^ are 

assumed to be 0 (m)-, at most (i.e., nearly quasi-steady operation), and 
solving for k, ^eld 

(ix. + JA) |=[| (1 + 2M) + Stl] ^ ^ m] + 0(m2) 

(16) 


a 


Thus, the frequencies of the possible standing mves are given by 

— = § (1 + 2N) + 6ti + 0 (m2) (17) 

The fundamental mode is ^ + 0(M^)^ which for quasi- steady 

operation becomes — = (The latter corresponds to the fundamental 

SL u 


^For nearly quasi-steady operation, is sm a ll and may arbitrar- 
ily be considered 0(M) even thou^ it is not functionally dependent on 
M. Therefore, if 0^^ = 0(M) , then = 0 (m) . If = 0(M), it then 

follows that = 0(M) (eq.. (9b)). 
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mode of an organ pipe, open at one end and closed at the other.) The 
system is unstable -when 

^ M + 0(M^) (18a) 

or 

^ =■ ^ M - SBi + 0 (m 2) (18b) 

■vdiich is consistent -with the original assumption that 
quasi- steady operation, the system becomes unstable ■when 

(^)i ^ ^ 

That is, the system is ims-table -when the slope of the characteristic 
curve exceeds the -value indicated by ecjuation (l9) . With increasing 
oiZi/a, the effect of 80^ must be considered. 

It may be concluded that the configuration defined by equations 
) and (l3c) ■will become unstable -vdien = 0(M), and the 

fundamen-tal frequency of the standing wave is — = 5 + 0 (m) . Increas- 

ing M or decreasing the slope of the inlet characteristic curve, dur- 
ing subcritical c^)eration, will tend to make a configuration more stable. 



%he characteristic curve of an inlet is frequently expressed as 
p^ against M or p-j. against W tdiere p^ is the singnatlon pres- 

sure in the combustion chamber and W = puS — is the corrected mass 

P«o 

flow referred to free-stream (or other) reference conditions. The rela- 
tion between and dp^/dM or ^ 


dW 


is given by 


f±.\ = (. 

[rvji du 


yp^ dM 


_1_ 


- M^|l + 0 (m2^ 


where it has been assumed that the s-tagnatlon temperature in the combus- 
tion chamber does not vary -('rith mass flow and that the ri^t-hand side 
of these equations are of order M (i.e., operation near point of insta- 

1 dp^ y dp^ _ 

= 0(m) and — = 0(M^) at the point of in- 


bility). Thus, ^ 
clpient instability. 


Pt 
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In appendix E, it is shown that the exit tends to danp incid.ent 
waves •vjhile the inlet anip3J.fies incident waves (for d±> O) • Thus 
the inlet is the component of the ram-jet configuration that is re- 
sponsible for instability. The role of the exit^ in dancing incident 
tra.ves, was previously pointed out in references 9 and 10. 


Burning Case 

The possibility of burning at station (sketch (a)) is now admit- 
ted. It is assumed that the configuration satisfies equations (?) and 
(l3)(with eq. (l3c) replaced by <<l). 

The boundary condition at Xg is, from equation (?), 


2 


K - 



( 20 ) 


while the botindaxy conditions across, the flame ffont are given by equa- 
tions (8 ) (with the heat release per unit mass assumed constant).' Take 
Xf to be the origin of the x-coordinate system (xj- = O), and define the 

foUoi/i°S constants: 


U = - 


U 


kxi 


B = 


- M^(l + - Mb(l +Mb)J 

With use of equations (6)> ?u,f/^u,f expressed as 

^u,f >, 1 - i cot U _iu' 

" 1 - i ^icot U " 


kxe 


B 


( 21 ) 


(22) 


It ^vas previously found that = 0 (m) and = 0(M) for the ram jet 
operating without biimlng at A => 0. It will be shown that Cj[ = O(Mj^) 
(and thus = O(Mu^)) for the coire spending burning ram jet. There- 
fore, to siniplify the present analysis, it is assumed, a priori, that 


^ = o(-%) 

Vf 


( 23 ) 
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\rfaich foUovra from equation (22), -vdien = O(Mjj), = 0(M^), and 

A » 0 (provided cot U << “). Using equation (23) changes the bound- 
ary Conditions across the flame front (eqs. (s)) to 





+ o(mJ) = 0 




^h,f 


(24a) 

(24h) 


Expressing aud Cg in terms of ^ and f ^<3. utilizing 

equations (20) and (24h) give 


t 


Mh 

^ (r^ - 1) + i cot B 


Mu 

1 + i (rt - i)cot B 


1 - 


"3 e 


IB^ 


-T (25) 


Tb - 1 C 2 cos B 


+ 0(M^) 


Substituting equation (25) into equation (24a) and expressing ^ ^ ic 
terms of give ’ 


Ci 


■(=v^) 




C^cot U + cot B b 


cot B 


C^cot B + cot U T 5 ■ 1 C]_cot B + cot U 

sin 


^3 

20 =:+^^ 


COS 


w 


5 C 2 cos B 


C^cot B cot U - 1 Mb C 3 b 

C^cot B + cot U 2 C 2 C^cot B + cot U. 




sin B 


+ 0(Mg) (26) 


Inhere 


r. 




'A 

1 

'Z 


In the derivation of equation (26) it was assijmed that 

1 C3 ei(B+B') 2 


Cb - cot U cot B 


- 1 c. 


cos B 


T ^5 


CbCot B + cot U 


< 0(1) (27) 
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Since the equation (27) is of order M for the nonhuming case, this 
assun^ition appears reasonable. Equation (26) indicates that 0^ = 0(M^)_, 

■which leads to = O(M^) and = O(M^) thus -verifies equation 

( 23 ) (provided that cot U « “) . Equation (26) is the cozmterpart of 
equation (l2). It defines "the acoustic iagoedance of a constant-area 
■tube ha'ving a through flowj a choked exit, and a planar flame front at 
X = 0. If an inlet is added at Xj^, -then k must have a value such 
•that the acoustic impedance of ■fche "tuhe matches the acoustic iiipedance 
of the inlet. Ihis condition .defines k and therefore -the s'fcahili'ty 
of ■the system can he determined. 

If it is assumed -that the flame is fixed (Av = O), and the gas 
constants are the same for both the burned- and ■uhbumed-gas regions 
(i.e., r^b = Tu^ ^b = ®u^ etc.), equation (26) becomes 




cot B + cot uV ^ 






cot U + 


. , ““(-A hJ cot b1 

5X-2 - (X - 1) ^ B J 


-^/r cot B cot 0-1- 




Bln B J y 


/■K)(l^) 

(E8) 


where X = equations such as ^ = -v/x, ^ = xjl + 0(M^)J , 

and ^ = -y/\ |l + 0(M^)J (based on eqs. (B2) of appendix B) have been 

used. Equation (28) has been solved numerically for -the following spe- 
cial case ; 


A = 0, ■»lr^ =» 0: Assume ■the configuration is operating at 1±ie 

neutral stability point (A = O), so that U and B are real. Fur- 
ther, assume that ■\}r = 0 (quasi-steady operation) . Equation’ (28) 

then becomes 


0 = -/X cot B cot U - 1 + 0(M^) 


(29a) 


«i 



(X- 1 ) 



(29b) 


Equation (29a) can be solved for o>. Using this -value of cD in equa- 
tion (29b) gives the ■value of 0^ at the neutral s-tability point. A 

ram-Jet configuration is "unstable ■when its -value of 0 jl is larger than 

that indicated by equation (29b). 
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Since U 
written as 


(i ^ 

■ au \ L / 


and B = — ^ (~ 


tA 


7i.,T \ Xg 

— 1 equation (29a) can be 
a^/ n 


toL /caL‘\ ^ ^ . _1 1 . 1 faiL\ ^ I 


+ 0(M^) 


(30) 


ojL 


which, for a given can be solved directly for — as a function of 

au 

This gives ~ as a function of ^ for the given X. Be- 
au 1* 

fine a frequency parameter 2 according to the relation 


Wj l; 


2 = = 2 


cdL 


yj 


(31) 


This parameter has the value of 1 for x^/L = 0 and Xg/L = l. (consid- 
ering the fundamental mode) and is somewhat larger than 1 for interme- 
diate values of Xq/L. For a given X, 2 has a single ma xi m um in the 


range 0 < -^ < 1 . 


This mftximnm can be shown to equal 
(■2)max = I tan-l(xlA) + 0 (Mu) 


(32) 


and occvirs at a value of Xg/L given by 


^ -.A 

^ + 1 


(33) 


(for X l) . Thus 1 <|j(2) mHv + 0(M.^2< 2 for K X A plot of 

2 against Xq/L is presented in figure 2 for X = 4 and 9. 


Equation (29b) is solved using the value of obtained from equa- 
tion (29a) . Substitution of equation (29a) into equation (29b) permits 
the latter to be written in the more convenient form 


% 


r-1 

2 


Mu 


1 + X cot^B 


lx + 


cos 


3X-2 - (X - 1) 


(vs O' 


cos B 


CO- 


t% 


+ 0(M^) 


(34) 
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For the nohbumirig case, 

( 34 ) indicates that — — 

2 


®i 




> 1 . 


1- When there is burning, equation 
Thus, a given configuration is made 


more stable because of the burning. Equation (34) gives 



- I/X. for 


Xq/l =a 0, and 



\ for Xq/L = 1. 


/r - 1 

The function ®i/'^~2 — ^ have many meixljtium and mlniTmTm points in 

the range 0 < Xq/L < 1 because of the relatively rapid oscillations 
of cos(B/-y\ Mjj) with Xq/L for Fana.ll A plot of g ^ 

against Xe/L for X = 9 and = 0-01 and 0.05 is presented in 


figure 3 . The value of ~ g " ^ increases from 1.89 at Xq/L = 0 

to 9 at Xq/L =1. Oscillations about the mean are more numerous for 
the = 0.01 case than for the = 0.05 case. 


These resiilts asstme a fixed planar flame front and a constant heat 
release per unit mass. The possibility of a combustion- induced insta- 
bility (AQ ^ 0 ) has not been considered. If any of these assumptions 
axe relaxed, considerably different stability criteria ■vroxald be expected. 


Ft 


CONCLUDIHG KEMAEKS 

A one -dimensional, analysis of ram-Jet buzz is presented. The par- 
ticular mechanism leading to btizz, which is treated herein, is the am- 
plification of aeoustic waves in the ram- jet combustion chamber due to 
successive reflections from the inlet and exit sections (i.e., lineax 
instability origin) . The system becomes unstable when the real part of 
the acoustic impedance of the inlet is greater than a term of the order 
of the combustion-chamber tfeich number. 

An analytical method for determining the acoustic impedance of an 
axbltrary inlet would be of interest. As mentioned previously, only the 
open- nose inlet has been treated analytically (ref. 9) and the analysis 
therein is restricted to low frequencies. 

It should be noted that the previous derivations are generally ap- 
plicable to configurations wherein the inlet of the present report is 
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replaced by a conipressor. It woiild appeax that only equation (9b) re- 
quires modification (to account for the variation of compressor dis- 
charge stagnation temperature with mass flow) . If the quasi-steady per- 
formance of the compressor is chareicterized by a polytropic exponent n 


[^1 = r4r “ M ,. -•) 


equation (9b) is re- 


placed by 


?1 



(35) 


for n of order 1, 1b at most of the order of and the solu- 

tions presented in the main body of the report (i-e.^ eqs. (15) to (34)) 
are applicable. If the flow in the combustion chamber can be considered 
isentropic, then appendix D app2d.es. 


Lewis F2JLght Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 30, 1955 
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APPENDIX A 


02,04 

°P 

Cv 

F,G,H 

k 

L 

I 

^±}T'q 

M 

N 

n 

P 

P(u) 

Q 

R 

S 


SYMBOLS 

The following symbols are used in this report: 
constant defining an^jlification 
speed of sound 

constants defined in eqs. ( 2 l) 
constants defined in eqs. (BS) 
specific heat at constant pressure 
specific heat at constant volume 
complex constants defined by egs. (4) 
functions defined by eqs. (3) 
complex constant (cd + iA) 
length of combustion chamber 
characteristic length 

characteristic length of inlet and exit section, 
respectively 

Mach number in combustion chamber 
integer ( 0 , 1 , 2 , . . .) 
polytropic exponent 

static pressure of gas in combustion chamber 

combustion-chamber static pressure as function of u for 
steady operation 

heat release per unit mass at flame front 
gas constant 

cross-sectional area of combustion chamber 
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T 

t 

U,U> 

u 

V 

W 

X 

r 

A 

5 

C 

T1 

e 

\ 

K 

p 

If 

a 

CD 

Subscripts : 


temperature of gas in combustion chamber 
time 

constants defined in eqs. (2l) 
axial velocity in combustion cdiafiaber 
velocity of planar flame ffont 
coinrected mass flow 

coordinate distance in stream direction 
ratio of specific heats 
net perturbation in a quantity 
perturbation frcm quasi-steady value 
acoustic impedance, 0 - iijf (eq. (5a)) 
entropy 

acoustic resistance 

ratio of tenq)eratin:es of-biimed to unbumed gases (!]}^J/!I^^) 


entropy impedance, — 

Cp ZjU 


(5b)) 


density of gas in combustion (chamber 
acoiistic reactance 

frequency parameter defined by eq. (3l) 
frecjuency, radians/sec 


b bumed-gas region 

e combiistion- (Chamber exit section 

f flame front 


i 


ccmbustion-chamber inlet section 
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r reference section 

u xmbumed-gas region 


I 
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APPENDIX B 

BOUNDARY CONDITIONS ACROSS PLANAR FLAME FRONT 

Consider a flame to be propagating with a small constant velocity 
V relative to a wall^ as indicated in sketch (b) . Conditions in the 


1 

^ 1 


i 

• 


7777777777777777 


(b) Flame moving relative to wall, 
burned- and unbumed-gas regions are related by 
Momentum; 

Pu + = Pb + pbK - 

State: 

Pu Pb 


= 1 


Continuity: 


PMi Pb^b^b 

PuK - = pbK - 


Energy: 




> (Bl) 


^b - Pb 
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■wiiere Q 1 b "tlie energy added, per unli; mass - If 'fclie flow velocities 
are gTnaii corapared with the speed of sound, equations (B1) become, 
respectively , 

^ = I + 0(Mg) 


■^u 

PuWi 


= I + O(I^) 


TuRuK : l) ^ ^ 

TbRbC^ - - 

where % = Uu/a^. Wow, assume that the flow is peiinirbed. The result- 
ing relations across the flame front are 


(B2) 




^ -^u 

^ ^b 




A(ub - v) 






(B3) 


From equations (id) and (2) _ 

^ = 1 
P r \V 

^ ^ 1 pa + (^ _ 1) 

T tL^ ^ Pj 

Substitution of equations (B4) into (BS) yields 


(B4) 


(B5a) 
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Auv 


Ut. 






,b 


^b% Pu/L 


0(1^)] 


(B5b) 


°p,t>^t) 1 
°P,U?T1 "'ft 


Anb 

°v,-b 


(Tb - 1) 


f£b 

Pb 


;]- ■ fes - ^) f [- ““«>] 


(B5c) 


Consider the eguations (B5) as applied to a ram- Jet combustion cbam- 
ber. Under steady operating conditions^ the flame front Is stationary 
(l.e., V = O) and Is located at the flameholder (x = x^)^ vhlch Is as- 


sumed to be planar. If sinusoidal perturbations are assumed, and these 
are expressed In complex form, division of equations (B5b) and (B5c) by 
(B5a) gives 




0 


(B6a) 


■where 





, ^ AQ 

Q ^,f 


+ 0(Kg) = 0 


(B6b) 



^2 


*^v,b^b 
• Cv,uTu 



= Tul 






The subscript f Indicates that -fche quantl-ty Is to be e-valuated at x^ . 
If "the heat release per ■unit mass Is assumed cons^fcant. Aft = 0 In equa- 
tion (B6b). The choice for the proper -values of Av/Zai^^f In equation 

(B6a) depends on the effectiveness of -fche flameholder In constraining 
the motion of the flame. For low-frequency oscillations, Av/Aa^^f can 

probably be taken equal to zero (l.e., the flame Is considered fixed at 
Xp). For higher frequencies, Av/Aa^^f probably cannot be taken equal 

to zero (despite -fche presence of the flameholder) and. In fact, may be 
complex . 
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APPENDIX C 


BOUNDARY CONDITIONS AT CHOKED EXIT 

AssTJme that the exit section is operating quasi steadily (coZe/a l) 
and is choked. Under these conditions, the Mach number at Xg (sketch 
(c)) is constant at all times, its value depending on the ratio of the 




(c) Exit section. 

tube cross-sectional area to the exit throat area. Thus (/5M)g 
. Using a2 = yp/p and equation (2) then yields 


0 or 




(Cl) 


From the definitions of ^ and §, equation (Cl) can be written as 




r - 1 




M(r 


^ 1 - 


(C2) 


which relates ^ and K at the exit. An expression s imil a r to equa- 
tion (C2) is derived in reference 9. 


If (uZe/a is not small, the appropriate boundary conditions at 
the exit must be determined by methods analogous to that used in 
reference 12. 
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APPENDIX D 


u 


SOLUTIONS OF EQUATIONS (ll) AND (l2) FOR ISENTROPIC FLOW 

In the main "body of the report_, solutions of equations (ll) and 
( 12 ) are presented, 'vdilch are applicahle for ram- jet configurations 
satisfying equations (l3). If it is assumed that the entropy waves 
in the comhustidn chamber are negligible (i.e., isentroplc flow) then 
^ = 0, and equations (ll) and (12) can be solved without recourse to 
the restrictions of equations (l3). The solution is presented herein. 

If 4 = 0, equations (ll) and (12) become, respectively. 


-2i 


kL 


a(l-M2) 


1 + ^i 1 - ^ 

^7741;: 


(Dla) 


and 


r - 1 


M + i cot 


Ci = 


C(1 ^M^)) 


1 + i 




(Dlb) 


Solving equations (D1) for k yields 


((D + iA)L It 
a(l-M2)“ 


f l(l + 2N) +- tan"^ 

2r l_02_ijf2 4 


In 


(l + 0i)2 + tf A' 


r-1 


(l-0i)2 + tf A + 


The frequency of the standing waves is therefore 


(uL 

a 


■I 2Tjf. “1 „ 

= i (1 + 2N) + i tan-1 i ^ (1 _ ^2) 

2 JC 1 - - ijr| 


(D2) 


(D3) 


If ijf^ and 9j_ axe functions of o>, then equation (D3) must be consid- 
ered as providing an implicit, rather than an explicit, solution for cd. 
V/hen i{rj^ = 0, the frequencies are given by — = -^(l + 2N)(l-M^). 
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Increasing the Mach, numher "tends "to decrease "the frequencies. The sys- 
tem will hecome uns"tahle ■^dleh 



If = 0 (hut M is not necessarily small) , equation (D4) becomes 



M< 0^ 


< 



(D5) 


which agrees wi"th pre"vious results for quasi- steady operation. If 
M « 1 (hut T}Tj^ is not necessarily small), equation (D4) becomes 

M(1 + ilff)[l + 0 (m 2)] < 0^ < [l + 0 (m2)] (ds) 

Equation (D6) indicates that a nonzero value of ijr^ will tend to in- 
crease the stab ili ty of a configuration, hut that the effect is negli- 
gible when iir| « 1. Note that equations (D5) and (D6) gi"v^ an upper 

hound, as well as a lower hound, to the -value of 0f which lead to 
ins"tahility. In all cases, positive -values of 0j_ are required. 

Equations (D1) to (DS) are valid pro-vided that the exit section 
is short compared -with the -wavelength of the oscillations and provided 
"that entropy -waves are negligible . There is no restriction on the rel- 
ative length of "the combustion chamber as compared with the inlet. 
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APPENDIX E 


SUPPLEMENTARY DISCUSSION OF RAM-JET STABILITY 

Additional insist into the source of ram-Jet instability can be 
gained by considering the ratio of the magnitxides of the reflected to 
incident waves at the inlet and exit of the combustion chamber. For 
convenience, the nonbiiming case is considered. 

Q?he pressure and velocity perturbations in the cconbustion chamber 
are related by 

Ap+ = pa An+ 

= -pa Au“ 

where the superscripts (+) and (-) designate quantities associated with 
dov/nstream and upstream traveling waves, respectively. The net pertur- 
bation at a point is given by 

(E2) 


ihi = Au+ + An" 


(El) 


If sinusoidal waves are assumed, and the perturbations are expressed 
as the real part of con^jlex quantities, then — ^ and equations 
(El) and (E2) give 


^ - 1 


(E3) 


The ratio of the magnitude of the downstream to the upstream wave at 
any point is 




C + 1 



c - 1 


(E4) 


The ratio of the magnitude of the reflected to incident wave at the 
inlet is then 




^i +1 

V^"ji 


^i-1 


(E5) 
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The ratio of the magnitude of the reflected to the incident wave at the 
exit is 



(E6) 


The condition that the system he unstable is that successive reflections 
of a wave (at the inlet and exit) should result in an increase in the 
magnitude of the wave. That is, the condition for instability is given 

t>y 



or, from equations (E5) and (E6) 


^i t 1 ^e - 1 
^i - 1 + 1 


> 1 


(E8) 


Since ^ = Q - i-^, equation (ES) can also be written as 


(0jL + 1)^ + tf (Oq - 1)2 + ^ 

(0^ - 1)2 + (eg + 1)2 + ijf2 


(E9) 


which agrees with the stability criteria indicated by equations (16) 
and (D2). 


For a short choked exit, 
so that 



2 

= ly Z . ~ i)m (iieglecting entropy waves), 
= 1 - (r - 1)M + 0 (m 2) (ElO) 


Thus, the reflected wave is weaker than the incident wave, and the exit 
tends to dan^p the wave system. Ihe larger the value of M, the greater 
the dancing effect. The system becomes unstable -vdien the inlet ampli- 
fies incident waves by a greater amount than they are damped at the 
exit. The values of 0^, for which this occurs, is given by equation 

(D4) . As indicated in appendix D, 0^ must be positive for the system 

to be -unstable. Thus, it is the inlet section that is responsible for 
the instability of ram- Jet engines. 


t- 

r- 

to 

to 
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(b) Subcritical operation. 
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(c) Typical perfomance curve. 

Figure 1. - Ram- Jet configuration and typical performance curve. 
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Figure 2. - Frequency pa 
(Eqs. (29a) and (31).) 
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